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We study the dynamics of a charge qubit, consisting of a single electron in a double well potential, 
coupled to a point-contact (PC) electrometer using the quantum trajectories formalism. In contrast 
with previous work, our analysis is valid for arbitrary source-drain bias across the PC, but is re- 
stricted to the sub-Zeno limit. We find that the dynamics is strongly affected by inelastic tunnelling 
processes in the PC. These processes reduce the efficiency of the PC as a qubit readout device, and 
induce relaxation even when the source-drain bias is zero. We show that the sub-Zeno dynamics are 
divided into two regimes; low- and high-bias in which the PC current and current power spectra 
show markedly different behaviour. To further illustrate the division between the regimes and the 
inefficiency of the detector, we present simulated quantum trajectories of the conditional qubit and 
detector dynamics. We also describe how single shot measurements in an arbitrary basis may be 
achieved in the sub-Zeno regime. 

PACS numbers: 73.63.Kv 85.35.Be, 03.65. Ta, 03.67.Lx 

I. INTRODUCTION 

Single shot quantum measurement of mesoscopic systems is recognized as an important goal. Fundamentally, it will 
allow us to make time-resolved observations of quantum mechanical effects in such systems, and practically it will be 
a necessary component in the construction of solid-state quantum information processors (QIP). There are numerous 
proposals for implementing QIP in solid state systems, using doped silicon^, electrostatically defined quantum dots^ 
and superconducting boxes^. In these proposals, the output of the QIP is determined by measuring the position 
of a single electron or Cooper pair. Single electrons hopping onto single quantum dots have been observed on a 
microsecond time scale using single-electron transistors (SET)^. Ensemble measurements of a double well system 
(qubit) have been demonstrated in superconducting devices^. So far, single shot qubit measurements remain elusive. 

It is therefore important to consider the measurement of single electron qubits by sensitive electrometers. Two 
possible electrometers have been discussed to date: SETs (see e.g. Makhlin et al.'^) and point contacts (PCs)^"^''. 
PCs have been shown to be sensitive charge detectors^^""'^^, and are the focus of this work. Figure 1 illustrates the 
physical system wc consider here, with a PC (represented by two Fermi seas separated by a tunnel barrier) in close 
proximity to one of the double well minima. Despite the apparent simplicity of this system it exhibits a rich range of 
behaviour. 

There are several energy scales of relevance: the splitting of the qubit eigenstates, 0, the bias voltage applied across 
the PC, eV = /xg — fin, and the measurement induced dephasing rate, Td, due to the different currents through the 
PC that result for different localized qubit electron states. These three energies define three distinct measurement 
regimes 

1. low-bias regime, where Td/'i ^ eV < <j>. 

2. high-bias regime, where Td/2 <C < eV. 

3. quantum Zeno limit, where (j) <C rd/2 ^ eV. 

In the quantum Zeno limit, frequent weak measurements localize the qubit, suppressing its dynamics^' To date, 
little attention has been paid to the low bias rcgime^^. Recently, the asymmetric power spectrum was calculated for 
arbitrary eV^'^ . This may correspond to the emission and absorption power spectra of quanta in the PC, as discussed 
in-^^. Unconditional calculations of the dynamics for arbitrary eV have also been discussed recently^^. 

We present a detailed analysis of the sub-Zeno limit, for arbitrary eV , showing the sharp difference in behaviour 
between the low- and high-bias regimes. Our analysis demonstrates the importance of inelastic tunnelling processes in 
the PC, in which electrons tunnelling through the PC barrier exchange energy with the qubit. These inelastic processes, 
illustrated in Fig. 1, incoherently excite or relax the qubit, yet are not in general experimentally distinguishable from 
elastic tunnelling processes using current measurement. 

The paper begins with a derivation of an unconditional master equation for the joint system consisting of qubit 
and PC, highlighting the significance of the Fourier decomposition of the interaction Hamiltonian. We then present 
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FIG. 1; Schematic of the qubit and PC showing lead energy bands. Electrons tunnelling from the source to the drain may do 
so elastically or inelastically, depicted by arrows. Different transitions induce different jumps, Pi, on the qubit. 

analytic solutions to this master equation, in both the high- and low-bias regimes. In Sec. V, we derive conditional 
master equations, which describe the dynamics of the qubit conditioned on particular realizations of the PC detector 
output. In Sec. VI we apply the results of the preceding sections to determine steady state properties such as the 
steady state density matrix (which demonstrates the close analogy between the PC bias voltage and a heat bath), and 
the corresponding detector output current. Following this, we derive steady state power spectra. We end the analysis 
with some sample trajectories, found by solving the conditional master equations, showing possible measurement 
outcomes. We conclude the paper with a discussion of the implications of our results for qubit spectroscopy, and for 
quantum information processing. 

II. SYSTEM HAMILTONIAN 

We model the double well system as a two level system, on the basis that the two lowest energy eigenstates, |e) and 
\g) of the double well system are well separated from higher lying single particle energy eigenstates. We write these 
eigenstates as 

|e) = cos(6'/2)|r) - sin(6'/2)|0, 
\g) = sin(0/2)|r)+cos(0/2)|O, 

where |Z) and |r) are the left and right localized states respectively. The Hamiltonian for the double well system, its 
interaction with the PC leads and the lead Hamiltonian is 

Hsys = -^a,-|a, = -<^a(^)/2, (1) 

^moas = ^{Tk^q+Xk,q<^z)ao,q^S.k + ^-C- (2) 

k,q 

Hicads = Hs + Hd = '^ujkial^k^s,k + a''^i^aD,k) (3) 
fe 

where = \l){r\ + \r){l\, a, ^ \l){l\ - \r){r\, = tan-i(f ), (j> = and tri"^ = \g)l^g\ - |e)(e|. We adopt the 

convention that Xk,q < so that the left well is nearest the PC. When = the energy eigenstates coincide with the 
localized states, since tunnelling is effectively switched off. At 6* = 7r/2, the tunnelling rate, A, dominates the bias, e, 
so the eigenstates are the completely delocalized states. 

III. UNCONDITIONAL MASTER EQUATION 

The von Neumann equation for the density matrix, i?, of the bath and system is 

fl(t) = -z[ifTot,i?]. (4) 

To derive the master equation for the reduced density operator, p, of the double well system we transform to an 
interaction picture with respect to the free Hamiltonian ~ Hgys + -ffioads, Hi{t) ~ e*^"*77Tote"*^"' — Hq = 
e'^^'i/nioasC"*^"*, expand the von Neumann equation to second order and trace over the lead modes, i.e. 

Pi{t) = TYs,D{-^[Hi{t),R{Q)] - f dt'[Hi{t), [Hi{t'),R{t')]]}, (5) 

Jo 
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where 



k.q k,q 

, f cos(0) e-'^'f' sin(9)\ 



3**"^ sin(6') - cos(6i) 

and Sk,q(t) is a time dependent system operator that may be written as a discrete Fourier decomposition Sk.q{t) ~ 
e~*^'^'=^'^''+'^"^*P„ for some time independent operators P„ and frequencies tj„. We can write the operator Sk,q{t) 
exphcitly as 

Sk^q{t) ^ e-i*('^'^-",)(e-"'^Pi + e"'^P2 + P3), (8) 

where Pi = ~ Xoo sin(^)|(7) (e|, P3 — (Pqo + Xoo cos(0)cri'^'') and we have assumed that the tunnelhng matrix 
elements Tk^q = Too and Xk,q = Xoo are constant over the regime of interest. 

The form of Sk,q{t) indicates that there are three possible jump processes, indicated in Fig. 1. P3 is associated with 
elastic tunnelling of electrons through the PC. Pi and P2 are associated respectively with inelastic excitation and 
relaxation of electrons tunnelling through the PC accompanied by an energy transfer (j). This energy is provided by 
the qubit which relaxes or excites in response. Inelastic transitions in similar systems have been described in^-'^, which 
calculated the current power spectrum through an open double well system due to shot noise through a nearby PC. 

Assuming the leads are always near thermal equilibrium, Tr s,D{o-i kP-i^)} " Tr5.£i{ai,A:P(i)} = and 
Tr5.D{a| ^,a„.fcP(t)} = Sijfi{uJk)p{t), where i,j G {L, R} and fi is the Fermi distribution for lead i. We also make the 
common assumption^'^ that if the lead correlation time is much shorter than other time scales, then the lower limit 
on the time integral in Eq. (5) may be set to —00. Equation (5) becomes 

pi{t) ^ - / dt' dujk / du;qgsiuJk)9Dit^q)fsi^^k){i ~ foi^^q)) 



— 00 



X {Sit)^Sit')pi{t') - Sit)pi{t')S{t'y - S{t')pi{t')S{t)^ + pi{t')Sit')^S{t)) 

- J dt' J dujk J dujqgs{i^k)gD{'^q)fD{i^q){l- fs{<^k)) 

X {S{t)S{t')^ pi{t') - S{t)^ pi{t')S{t') - S{t')^ pi{t')S{t) + pi{t')S{t')S{t)^) , (9) 

where we have made the standard replacement — > J dLUkgs{D)i^k) where gs(D){^k) is the density of states for 
lead S{D) which we will hereafter assume is a constant, gs(D)i iii the energy range of interest. 

We wish to evaluate the integrals in Eq. (9) at zero temperature, so to show the method and approximations used, 
we evaluate a particular term in Eq. (9) using the Fourier decomposition of S{t): 

J2fs{cuk){l-fD{iOq)) f dt'S{t)pi{t')S\t') 

k,q 

= Y,U^kl diUggsgDfs{^k){l-fD{iOq)) f dt' e-'(^--^-+^-^'e'^^--'^-+^-^'' P^pi{t')Pl 
~ ^gsgoY, /d^fc fdujgfsMil - fD{uJq))5{^k - +t^„)e*(""-"'")*P™pi(t)P,t, 

~ T^gsgO ^ / dcjfc/5([Jfc)(l - foii^k + W„))P„/9l(t)P|. 

ri 

= T^gsgO J2 - + ^n)PnPlit)Pl (10) 

n 

where Q{x) = (x+ |a;|)/2 is the ramp function and pi is the chemical potential of lead i. The first equality follows from 
substituting the Fourier decomposition of S{t), the second equality follows from making the Markov approximation, 
pi{t') — > the third (approximate) equality we have made a rotating- wave approximation (RWA), where we make 

the replacement g^{'^r^~^m)t _^ (5,„ „, and finally we integrate over iOk- The RWA is reasonable as long as ^ v^eV, 
where ly = y^'^T^gsgoXoo is a small quantity. This is justified when we come to integrate the master equation where 
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terms with n ^ m are rotating sufficiently rapidly to vanish. Using the same arguments on each term in Eq. (9) 
results in the general form of the master equation 

= 2TigsgD {vyQ{eV + uj^)PMt) + Q{-eV - u,,)Pl]p,{t)) , (11) 

n 

where 1^ = /is — hd is the bias applied across the PC, 'D[B]p ~ J[B]p — A[B]p, J\B]p = BpB^ and ^[-B]p = 
^{B^Bp + pB^B). 

We may transform the master equation, Eq. (11), back to a Schrodinger picture to reinstate the free dynamics 
of the double well system. It conveniently turns out that the transformation of the P„'s back to the Schrodinger 
picture merely multiplies them by the unitary factor e*""*. The Lindblad superoperators are invariant under such a 
transformation, so in the Schrodinger picture 

m = -^[i/sys, p{t)] + J2 + L^nMcnW) + H(-ey - Lu,,)V[ci]p{t)) = Cp{t), (12) 



where H(x) is the unit step function. For convenience, we have defined c„ = -^Z Q{eV + uJn)Pn, so that 

ci ^ ,y^eVT^sm{e)\g){e\, (13) 

C2 = iy^\eV-cl,\sm{0)\e){gl (14) 

C3 ^ iyVeVcos{e)ai''^ +TVeV. (15) 



We will show later, using a explicit model of the detection process, that T = y/2ngsgjjToo . 

IV. ANALYTIC SOLUTIONS OF THE UNCONDITIONAL MASTER EQUATION 

A. High bias regime 

In the high bias regime, H(— eV^ — ujn) — for all ri, and the interaction picture master equation takes the form 

pi{t) =Y,T^[c,,]pi{t). (16) 

n 

We can solve the unconditional master equation in this regime exactly. We present the solution in the interaction 
picture, since it contains the relevant decay time scales. (It is straightforward to find the corresponding dynamics in 
the Schrodinger picture via the transformation p{t) = e~'^=»'=*pi(f)e*^"y=*). The general solution is given by 

p^,(t) = e-r^(i+-^^W)*/2p,,(0), 

- ^-e— ^^^^*(^-p..(0)), (17) 

where we have defined Td = 2i/^eV, consistent with^. 

B. Low bias regime 

Our analysis is also valid in the low bias regime, (j) > c-V. The unconditional master equation in this case, in the 
interaction picture, may be written 

piit) = V[ci]piit) + V[cl]piit) + V[c3]piit). (18) 

Again, we can solve the unconditional master equation in this regime exactly, and we give the solution in the interaction 
picture, 

P^t) - e-(i^--^(«)+^^'^^'"^W)*Pse(0), 

p,,{t) = l-e-2''^^-"^W*(l-p,,(0)). (19) 
Note that at (j) — eV, Eqs. (17) and (19) agree. 
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V. CONDITIONAL MASTER EQUATIONS 

It is also possible to describe the dynamics of the continuously measured charge qubit using a conditional master 
equation (CME). This describes the evolution of the state of the system, conditioned on the previous detector output. 
The CME is useful, since it allows individual measurement runs to be simulated, and also allows a straightforward 
calculation of detector output statistics, such as the average current and the power spectrum of current fluctuations. In 
this section, we derive a CME using an explicit model of the measurement process in terms of projective measurements 
of the number of electrons that have tunnelled through the PC. Our derivation follows a similar argument as that 
presented in^'* for the case of a continuously observed fluorescing atom. 

Owing to the different types of allowed jump processes in each regime, we treat the high and low bias regimes 
separately. 

A. High bias regime 

In order to derive the CME in the high bias (eV > 4>) regime, we first consider the unconditional evolution of the 
qubit over a short time interval (5i, which is taken to be sufficiently short that no more than one electron tunnels 
across the PC during this interval. This unconditional evolution can be written as a sum of two conditional terms, as 

Pi{to + St) = poc{to + St) + Pic{ta + 5t) . (20) 

Here, poc is the unnormalized conditional density matrix for the qubit, corresponding to having observed zero electrons 
pass through the PC in the interval 5t, whereas pic is the unnormalized conditional density matrix for the qubit, 
corresponding to having observed a single electron pass through the PC, in the source to drain direction, in the same 
interval. In the high bias regime, tunnelling in the reverse direction (drain to source) does not conserve energy and 
therefore is strongly suppressed on timescales of interest, as we discuss further below. The normalization of these 
conditional states is chosen such that Tr {po(i)c(^o + St)^ is the probability that zero (one) electrons tunnelled through 
the PC in the interval St. 

Assuming that at time t^, the total density matrix for the system and detector may be written in the factorized form 
p{tQ) ® |0)(0|, where |0)(0| represents the zero temperature state of the source and drain leads before any electrons 
have tunnelled, pic(io + St) may be written 

Pic(io + St) = Trs,D {niUiito, to + St) {p{to) ® |0)(0|) uj{to, to + St)Ui^ . (21) 

Here, 

E ^D.a'^sJO){0\4,(3<^D,o. (22) 

is the projector onto the subspace where only one electron has tunnelled through the PC in the source to drain 
direction, Ui{to,to + St) is the interaction picture time evolution operator for the interval St and fc^j^j) is the Fermi 
wavenumber for the source (drain) lead. For brevity we will suppress the range of the summation below. Note that 
Eq. (21) follows directly from the projection postulate of quantum mechanics. 

An approximate expression for pidto + St) can be found by first expanding the time evolution operator up to second 
order in St, i.e. 

Uiito,to + St)Kl-i dt'Hi{t')^ j dt' dt"Hi{t')Hi{t") , (23) 

J Iq J to J to 

and then substituting into Eq. (21) to give 

Picito+St) = J^i^l'^kpao^aUiita^to + St)\0)p{ta){0\ujito,to + St)a^^^^asM, (24) 

a, 13 

J2 / dt"{0\al^aD,^Hiit')\0)pito){0\Hiit")a^^^^as^p\0). (25) 

a,/3 " " 

Note that the zeroth and second order terms in the expansion of Ui{to,to + St) do not contribute to the above 
expression for pic{to + St), since the expression contains a projector onto the subspace where only one electron has 
tunnelled through the PC. 
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Evaluating the integrals in Eq. (25) leads to the expression 

Plcito +St)=J2 [fo^PlPl + + fo^M p{to) [UfnPl + U132P2 + , (26) 

where Pi are the operators introduced in Sec. Ill, and 

/a/31 = St e-'(-^— ^+'^)(*o+^*/2)sinc ({uj(3 - ^a + cb)5t/2) Hifis - Lo^)H{u,3 -pd), (27) 

/a/32 = 5te-'('^''-"°-'^'(*°+^*/2)sinc((t^0-^„-</))5i/2)if(Ms-^a)i/(c^/3-AiD), (28) 

/a/33 = St e-*(-^— »)(*o+*t/2)sinc {{ojp - Lo^)St/2) Hifis - Lo^)H{ujp - fio) ■ (29) 

For values St ^ (f)^^ , the cross terms in Eq. (26) make a very small contribution, and can be neglected. The 
remaining terms in Eq. (26) may be evaluated by converting the sums into integrals via the replacement — s- 
/ d(jJkgs{D){^k) as in Sec. III. Finally, again for St ^ 4)^^ , wc make the following replacement in the integrands 

sm^{xSt/2)/x^ 7rStS{x)/4:. (30) 

Thus Eq. (26) reduces to 

3 

Plcito + St) ^StY,JK]p{to), (31) 

n=l 

where the jump operators {c„} are given by Eqs. (13 - 15), with T = y/2TFgsgDToo. 

Note that using the replacement of Eq. (30) amounts to a statement of energy conservation (on sufficiently long 
timescales). For this reason, processes such as reverse (drain to source) tunnelling are suppressed in the high bias 
regime, since they do not conserve energy. 

The unnormalized state conditioned on having observed no tunnelling events in the time interval St, poc, can be 
found by repeating the above derivation, but replacing Hi in Eq. (21) with IIq = |0)(0|. However, it is also possible 
to find pqc by a more direct route, by noting that the normalized, unconditional state satisfies 

piito + 5t) =(^1 + StY,T^[cn]j Piito) ■ (32) 

Combining Eqs. (20), (31) and (32), gives 

Pocito + St)= (^l - St Y,Acn]j Piito) ■ (33) 

Equations (31) and (33) can be combined into a single, normalized, CME for the conditional state of the system, 
in the Schrodingcr picture, as 

dpcit) = dNit) _ ij p^{t)+dt (^-Y,^A[c,,]pcit) + Picit)pcit) - i[H,ys,Pcit)]j ■ (34) 

Here, pdt) denotes the normalized conditional state of the system at time t, and Pidt) = Tr |X]ri=i »^[cn]pc(i)| is 
the probability density of observing a single electron tunnelling event. We have also introduced the classical stochastic 
increment dNit) € {0, 1} which denotes the number of electron tunnelling events in the interval dt. The expectation 
value of this increment is given by 

E[dNit)] = Picit)dt = Tr J[c„]pe(t)| dt . (35) 

Note that although the CME of Eq. (34) has been written in terms of an infinitesimal time increment dt, it should 
be understood that in deriving the CME we have made use of a coarse- graining in time. In particular, neglecting the 
cross terms in Eq. (26) means that the derived CME neglects dynamics on fast timescales of order 0^^. Implications 
of this will be discussed further in Sec. VII C. 
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B. Low bias regime 

111 the low bias {eV < (f>) regime, the CME must also take into account tunnelling in the reverse (drain to source) 
direction. Note that these reverse tunnelling events are, in principle, distinguishable from the forward tunnelling 
events. Thus the unconditional evolution over the interval St may be written as the sum of three conditional terms: 

Pi{to + St) = pocito + St) + pi+c(to + St) + pi-rik) + St) , (36) 

where pi+c (pi-c) denotes the unnormalizcd conditional state of the qubit corresponding to having observed a single 
tunnelling event in the forward (reverse) direction in the interval St, and poc corresponds to having observed no 
tunnelling events in the same interval. 

The derivation of the CME now proceeds in the same manner as in the high bias case, except that now three 
projectors are used, corresponding to the different distinguishable measurement outcomes: 

Uo - |0)(0|, (37) 
ni+ = ^a\y,^as,p\Q){Q\alfjaD.a, (38) 

a, (3 

= E4,a«^,/3|0)(0|aL,/3as,a. (39) 

a, (3 

The resulting CME may be written 

^ {"^^ - o-W + '"v-«) - i) Mt) 

+ dt (- {.4[ci] +yl|4] + Alcl} pAt) + [PuJt) + Pi-M Pcit) - ![//.,,, . (40) 

Here, dN^{t) G {0,1} denotes the number of forward (reverse) electron tunnelling events in the interval dt. The 
expectation values of these increments are given by 

E[dN+{t)] = P,+,{t)dt = TT{J[cMt) + J[c3]Pcit)}dt, (41) 
E[dN~{t)] = Pi-c(t)fit = Tr|j^[4]pc(t)}rft- (42) 



VI. STEADY STATE RESULTS 

A. Pseudo-thermal ground state occupation 

We can use the unconditional master equation in each regime to compute the unconditional, steady-state probability 
of the qubit to be found in its ground state. It is given by pgg(oo) = {g\p{oo)\g). From Eqs. (17) and (19), we see 
that pgg{oo) is 

f ^ / ^ if</)<ey 
P..M = | Y iil>eV 

In contrast, the ground state occupation of a qubit in thermal equilibrium with a heat bath at temperature T is 
given by p*'^°''™(cx)) ~ j^'^^pr- Figure 2 shows these two ground state probabilities for the case where ksT — eV/2, 
and there is an evident analogy between the PC bias voltage and an external heat bath, even though the leads are 
each nominally at zero temperature, that is /s,d(w) = II(/xs^d — uj). As the PC bias gets larger the relative fraction 
of each component of the mixture tends to 1/2, analogously to a qubit in contact with a heat bath at very high 
temperature. 

It has been suggested elsewhere that the PC bias acts somewhat like a heat bath^^, and we note that this corre- 
spondence is a natural conclusion of our analysis. It may be seen that even in the limit eV 0, the PC still causes 
the qubit to decay to its ground state. We can see from Eq. (19) that in this limit, the energy relaxation time is 
given by = 2i^^0sin^(0), and the dephasing time is rj"^ = r(jCos^(6') + ;^^0sin'^(6'). When the PC bias voltage is 
zero, we see that = 2t,-, indicating that the dephasing is solely due to energy relaxation. This zero-bias dephasing 
arising from relaxation is the solid state analogue of the optical decay of an atom in a vacuum. 
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FIG. 2: Equilibrium ground state occupation probability for a qubit near a PC (solid), and the thermal equilibrium ground 
state occupation probability, pgg'^''™ (oo) for a qubit in contact with a heat bath at temperature T = eV 12 (dashed). 



This is a significant practical issue, since it means that the PC cannot be turned off merely by making the PC bias 
zero. To perform coherent quantum logic operations on the qubit in the presence of the PC, with error rates small 
enough to permit fault tolerant quantum computation, will require an extra gate to control the coupling, between 
the qubit and the PC. v can be significantly reduced by raising the height of the tunnel barrier between the source 
and drain leads of the PC. Then when the PC is off (i.e. eV = 0), v'^(\) can be made much smaller than by a factor 
of order the fault tolerant threshold^^ ~ 10~^, by raising this barrier. 



B. Steady state current 



The steady state current can be calculated using expectation values of the stochastic increments defined in Sec. V. 
In the high bias regime, the steady state current is given by /ss = eE\d'N [t) j di\t-,a^ = eTT{J2n=i '^[^n]p{<x)}. In 
the low bias regime, it is given by Iss = eE[{dN^{t) — dN^ {t))/dt]t-^oo = e Tr{ J'[ci]p(oo) + J'[c3]/9(oo) — J[cl]p{oo)}. 

The current can be conveniently expressed in terms of the expected currents, // and 7^, for the case when the 
tunnelling between dots is switched off {9 = 0), and the qubit is in the localized state |^)(Z| or |r)(r|, respectively. 
Then the average value of these currents is given by Jioc = ^^-^^ = e(T^ + iy^)eV and the difference is (5/ioc = I,- — h = 
-UvT eV. 



To first order in (5/ioc, we find 




cos{e)4> gj^^^ ^ j^^^ _ e jjj^^ if < 



2eV ""^ioc — -Moc 2eV' 
^^^SIioc = Iloc - J^SIioc if (/)> eV 



(44) 



Using Iss we can compute the DC conductance 

( T^ + ,.^ + 2iyTe/{eV) ii ^ < eV, , 
^-^''^^ ~ \ T^ + iy^ + 2i^Te/(b if0>el/, ^^^^ 

in units of e'^/h. Note that these results agree with those found in Rcf. 26, and given implicitly in the shot noise 
component of the power spectra of Ref. 19. 

Equation (44) indicates that the steady state current depends on the steady state charge distribution of the qubit, 
which in turn depends on the qubit bias e. As e varies the conductance changes accordingly, which is qualitatively in 
agreement with observations of charge delocalization made in recent experiments^^'^^'^^'^*. This behaviour is shown 
in Fig. 3. Note that the cusp at e = a/ (eV)"^ — appears when eV > A. 

This result suggests an improved method for performing spectroscopy of the qubit using steady state conductance 
measurements through the PC: the value of the applied voltage, eV^ at which a cusp in the conductance (when plotted 
as a function of e) becomes evident directly gives the unbiased qubit tunnelling rate, A. Alternatively, observing the 
conductance as a function of eV, for a fixed value of e, yields the splitting 0, by determining the voltage at which 
the conductance changes from low- to high-bias behaviour in Eq. (45). These methods offer a potential improvement 
over existing, indirect techniques^*, in which an estimate of A is extracted from a fit to a model, which in turn is 
limited by the uncertainty in the estimated electron temperature. 



FIG. 3: The steady state conductance, G, versus t/eV , for different values of the double well tunnelling rate, (light) eV/ A = 0, 
(medium) eV/A = 2 and (dark) eV/A = 4 



VII. POWER SPECTRA OF CURRENT CORRELATIONS 

Using the conditional master equations, we are able to evaluate the steady state power spectrum of the two time 
current correlation function 

G(t) =E[I{t + T)m] -E[I{t + T)]E[m] , (46) 

and we adopt the same procedure as described in the appendix of Goan and Milburn^. The power spectrum is the 
Fourier transform of this quantity, 

dTG{T)e-''^^ = I dTG{T)cos{ujT), (47) 

-OO 

since G(t) = G(-t). 







A. High bias power spectrum 

To compute the steady state correlation function, we take the limit t oo, so p{t) poc- Using the relation 
I{t) dt = e dN{t), where dN{t) G {0, 1} is the number of electron tunnelling events in the time interval {t, t + dt) and, 

E[dNit + T)dNit)] - Proh[dN{t)^l]E[dN,it + T)UNit)=il 
Prob[diV(t) = 1] = dtTi-{pic{t + dt)} , 

E[dN,{t + T)UNit)=i] = dtTi-{J2j[cn]E[pic{t + T)Ur4it)=i]}, 

71 

E[pic{t + T)UN(t)=i] = e^^pic{t + dt)/Tr{pi,{t + dt)}, 

dN{tf = dN{t), (48) 

we obtain the correlation function in the high bias regime for t > dt 

e^Tr{j:^^,J[c^]e^-J[c^,]p^}~TT{j:^J[cM') for r > dt, 
^hb(rj - I ^2 ^^^^^ J[c„]po.}Sir) for r = dt. ^^^> 

We write Ghb(''') in terms of /loc and SIioc'- 

Ghb(r) = e(/ioc - WiocTT^ cos{e))S{T) + e~^-^'''"^''^-e'{l - .^y{2TeVcos{e) + z.0sin2(0))2. (50) 

Dropping the small term j/0sin^(0), which is formally equivalent to keeping only the lowest order term in a series 
expansion in i5/ioc, and taking the Fourier transform gives 

.... ,j , <^/Lr4i-(^)Bin^(2e) 

S,,[.) ^ So - eM..- cos(^) + 4(r2,i,4(,)^^,) (51) 

where Sq = 2e/ioc- 
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FIG. 4: Power spectra, S{ui)/So, computed in this paper using Eqs. (51) and (54), as a function of eV/<j} for S = 7r/4, i^^ = 0.02. 

B. Low bias power spectrum 

The correlation function in the low bias regime turns out to be time independent. In this regime, there are two 
distinguishable jump processes that can occur: source to drain electron tunnelling (both clastic and inelastic ), and 
drain to source tunnelling (inelastic only). Therefore the current is related to the number of jumps in the time interval 
by I{t) dt ~ e{dN^{t) — dN^{t)), where dN'^'-^'>{t) counts the number of source-to-drain (drain-to-source) tunnelling 
events in the time interval {t, t + dt). 

To evaluate the expectation values in G(r), we use slightly generalized versions of Eqs. (48), for instance 

PToh[dN^{t) ^ 1] = dtTT{pi±^{t + dt)} , 
E[dN^{t + T)UNy(t)=i] = dtTT{J-E[pic{t + T)UNy(^t)^,]}ioTx,y==±. 

As before, to compute the steady state power spectrum, we make the replacement p(t) p^o = l5)(5|j for the low 
bias regime. It is evident that Poo = J\c^-2\Pcx, = J[ci\Poc ~ 0, so we find that 



Gib(r) = 




= e^{r + vcofi{e)YeV5{T) (53) 

Thus, the steady state low bias current correlation function is just due to clastic tunnelling through the PC. The low 
bias power spectrum is 

5ib(a.) = 5o - e «ioc cos{e) + 0{5ll,). (54) 
The power spectra computed in this paper are shown in Fig. 4 for = 0.02. Note the sharp change at eV = (p. 

C. Comparison with previous results 

The work presented here is valid when <C (f), whilst previous work has dealt with the large bias limit, ^ eT^^'^'^*^. 
These limits are shown schematically in Fig. 5(a). The limit F^ <C ^ eV is common to both this work and previous 
analyses^'^'^'', so we can compare the predictions of each in this limit. 

The current correlation function, Gqm{t), given in^ (subscript GM denotes the authors initials) is only valid for 
e = 0, but it is straightforward to calculate the general result, which is useful in order to compare our results with 
previous work. We present the general calculation in Appendix A. Figure 5(b) shows the power spectra computed in 
this paper (dark) and using Ggm{t) (light) for large bias, ^/eF ~ 0.01, and weak coupling, F^ = 0.1. The central 
peak at w = 0, due to inelastic transitions between the qubit energy eigenstates, shows good agreement between the 
results. There is a marked difference between the spectra at frequencies near In particular the peaks predicted 
in^ are absent from our results. However, this difference can be understood by noting that, in deriving our CME in 
Sec. V, we have employed a temporal 'coarse graining', which implies that our CME does not describe dynamics on 
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(a) ■ (b) 

FIG. 5: (a) Phenomena predicted in this paper (sub-Zeno) and previous work (large bias), (b) The power spectrum computed 
in this paper with Vd = 0.1, 6 = 7r/4 and eV = 100(/f> (dark), and the power spectrum computed by®, (see appendix A), with 
the same parameters (hght). 

short timescales of order 4>^^ . Therefore, our expressions for the power spectra (which are derived from our CME) 
are vahd only for frequencies lo <^ (p. 

In the low bias case, we predict a completely flat power spectrum (Eq. (54)). Although this prediction is also based 
on a coarse-grained CME, there is reason to expect that this should be valid at all frequencies. In particular, since 
the qubit rapidly relaxes to the (pure) ground state \g), which is a stationary state, there should be no oscillatory 
signal in the PC current, and no peaks in S'ib(i^) should be seen ai uj ~ cj). This prediction is in agreement with those 
made independently in Ref.^". 

VIII. MEASUREMENT TIME 

The detector partially projects the qubit onto the energy eigenbasis, as seen in the form of jump operator C3. The 
energy eigenbasis is therefore the 'preferred' basis for the detector. 

An important quantity of interest is the measurement time. This is the time taken for the measurement to project 
a qubit initially prepared in an equal superposition of energy eigenstates onto one or other of them. We therefore take 
p(0) = (|e) + |g))((e| + (ff|)/2. Following Goanetal.^ we compute the rate of change of £;[dz2(0)] = E[TT:{a^^^ dp^{0)Y]. 

We present the calculation for high bias, and note that we obtain the same result for low bias. Using the fact that 

dpc{t) ^ pcit + dt) - pc{t) = dN{t)pic{t + dt) + (1 - dN{t))pocit + dt) - p^(t), (55) 
dN{tf = dN{t) and Tr{cri'''p(0)} = 0, we find that to first order in dt 

E[dzlm - i?[diV(0)Tr{a(^)pi,(dt)}2]/2, 

2Tr{pie(dt)} ' 

= TdCos^{e)dt + 0{v^<p). (56) 

Thus the measurement time is t^-^^ = TdCos^{9) + 0{v^(t)). We note that t^i > t^, where Td is the dephasing time 
for decay of the off diagonal elements in Eqs. (17) and (19), with equality only when = 0. This indicates that the 
detector is inefficient unless the qubit energy eigenstates are localized. 

IX. QUANTUM TRAJECTORIES 

In Figs. 6(a) and 6(b) we show two different quantum trajectory simulations for 6 — 0, using the same parameters 
to generate both (see caption for values). We show the conditional probability (dark, solid line), {l\pc{t)\l), and 
unconditional probability (dark, dashed line), {l\p{t)\l), of the qubit to be in the left well. Also shown is the purity of 
conditional state (hght, solid), Tr{pc(0^}: of the unconditional state Tr{p(t)^}, (light, dashed). We also show 
a histogram of the number of jump processes that were encountered in each time window, and this corresponds to 
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FIG. 6: (a) & (b) Showing distinct simulations where the qubit collapses to |/) and \r) respectively [for 9 — Q, eV/<f) — 10, — 
0.005, = 0.5]. At = 7r/2 currents are uncorrelated with the qubit in either (c) high-bias [eV/(l) = 10, = 0.005, = 0.5] 
or (d) low-bias [eV/4> = 0.5, i^^ = 0.05, — 0.5]. In all panels, dark curves are {l\pc{t)\l) (solid) and {l\p{t)\l) (dashed), and 
light curves are the conditional purity Tr{pc(i)^} (solid) and unconditional purity Tr{p(t)^} (dashed). Histograms show the 
number of jumps per time interval, (scaled by 1/10). 



experimentally measured currents. In all trajectories shown in this paper, the histograms are scaled by the factor 
jq, and are mutually comparable. In these two figures, the initial state of the system is chosen to be the pure state 
(IO + k»/^/2. 

The evolution of the measurement clearly shows the qubit collapsing to the state |^) or |r), and the state remains 
pure throughout. This is because, at 6* = 7r/2, the inelastic jump operators ci and C2 are suppressed, and the remaining 
jump operator C3 partially projects the state onto the localized basis. The average current, determined by the mean 
of the histogram is higher in the case where the qubit collapsed onto the state |r), which is consistent with physical 
intuition, since in this configuration, the tunnelling rate through the PC is highest. The fact that the unconditional 
probability to find the electron on the left well is constant at 0.5 reflects the fact that the state of the system collapses 
to \l) or \r) with equal probability, consistent with the initial state preparation being an equal superposition of the 
localized states, and thus for = 0, the PC serves as a good quantum non-destructive (QND) measurement device. 

Figures 6(d) and 6(c) shows the dynamics of the system for 6 = 7r/2 in the high bias {eV/4) = 2) and low bias 
{eV/cf) — 0.5) regimes, respectively. In both cases, the initial state of the system is |^). 

As predicted by Eq. (17), in the high bias limit the equilibrium state of the system is an unequal mixture of \g) 

and |e). We may calculate the unconditional steady state purity, and we find Tr{/9(oo)} = '''^'^2e\rf = 0.505 for the 
parameters used. This agrees with the simulation of Fig. 6(d), where the conditional and unconditional evolution is 
very similar. 

Similar comments apply to Fig. 6(c), however the unconditional steady state is always \g) for the low bias limit, which 
is of course a pure state, and for the case considered is just an equal superposition of localized states \g) = (\l) + \r))/y/2. 
Thus the probability for the qubit to be in the left well is just the same as for Fig. 6(d), but the equilibrium state is 
pure. 

We illustrate the effect of changing parameters in Fig. 7. We show simulated trajectories for five different bias 
voltages, eV = 0.5, 1, 2 and 4, and for three different mixing angles 9 — tt/8, tt/4 and Stt/S. In the high bias regime, 
and small values of 9 there are fluctuations that bear a superficial resemblance to the quantum Zeno effect (QZE). 
The state of the qubit tends to fluctuate between mixtures of energy eigenstates, and the conditional density matrix 
is diagonal in the energy eigenbasis, with off-diagonal terms being negligible in this basis. 

This resemblance to the QZE is only superficial however. In the QZE, the free dynamics of the qubit are suppressed 
due to the strong and frequent measurement which project the system into pure, localized states {|Z), |r)}. The sharp 
transitions that occur in the QZE result from the small probability for the electron to make a transition between sites 
which is realized as a rapid, occasional tunnelling event. 

In contrast, the transitions that are evident in Fig. 7 for small 9 are due to relaxation or excitation of the qubit, 
accompanied by inelastic tunnelling processes through the PC, rather than the 'collapse of the wavefunction' due to 
accumulated information about the qubit in the case of the QZE. As 9 increases, the transitions become smaller in 
amplitude and less distinct, until 9 = 7r/2 when inelastic tunnelling events through the PC dominate the current and 
the qubit state is damped to the completely mixed state. 

At the transition from low bias to high bias, eV/(p = 1, fiuctuations begin to appear in the steady state trajectory. 
For eV < (j) the conditional steady state of the qubit is constant, whereas for eV > cj) the conditional steady state 
shows fiuctuations which increase with increasing bias voltage. 
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e = 3n/S 




FIG. 7: Simulated trajectories from left to right, 9 = tt/8, 7r/4 and Stt/S. From top to bottom, eV/(j> = 0.5, 0.9, 1.1, 2 and 4. In 
all = 0.05 and = 0.5. In each panel, the dark curves are {l\pc{t)\l) (solid) and {l\p{t)\l) (dashed), and the light curves are 
the purity of the conditional state Tr{pc(t)^} (solid) and the unconditional state Tr{p(t)^} (dashed). The histogram shows 
the number of jumps that occurred in each time interval. The histograms are all scaled by the same factor, so are comparable 
with one another. Reverse tunnelling events are responsible for the sharp changes in the first two panels. 



X. DISCUSSION 



The predictions of this work should be experimentally verifiable. In particular, it should be possible, in near-future 
experiments, to observe the qualitative change in both the steady state current, and the fluctuations in the current, 
in the transition from the low- to high-bias regimes. This sharp transition between low- and high-bias regimes also 
provides possible techniques for performing spectroscopy of solid state qubits, as discussed in Sec. VI B. 

The results presented here describe the dynamics of the measured qubit in the sub-Zeno limit, Td ^ 0, for arbitrary 
PC bias eV, whereas previous analyses^'^'^"'^"^ are valid only for large PC bias, eV ^ cf>, but arbitrary Td/(j)- Outside 
these regimes, for finite ratios Tdl 4> 3'iid eV/(j) neither approach is formally valid, and non-Markovian effects may play 
a significant role. This parameter regime is therefore an open area for investigation. 

There are a number of practical issues that the present paper raises. Firstly, as discussed earlier, to turn the 
measurement off, it is not enough just to turn the PC bias to zero. Simultaneously one must also make small, 
which can be accomplished with extra surface gates. 

Secondly, in order to perform a good single-shot measurement we require that the qubit measurement time is 
much shorter than the detector-induced relaxation time,'^ i.e. r„i <^ t^. From Eqs. (17), (19) and (56) we see that 
— 2iP max {(j),eV)sm^{e) and t'^ = 2TdC0S^{e), so we require that sin(6') w 0. Therefore, our work indicates 
that measurements in the sub-Zeno regime are not possible for highly delocalized qubit eigenstates. Furthermore, 
ideal measurements (i.e. measurements for which the conditional state remains pure) are only attainable when 6* = 0, 
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so the energy eigenstates are the locahzed states. 

In order to perform good single-shot measurements, one could operate in the Zeno regime, though there may be 
technical difhculties in obtaining a sufficiently large value of r^. since this requires passing large currents through the 
PC. These large currents can, in turn, lead to heating in the detector, which can significantly increase noise and lead 
to additional decoherence of the qubit.^^'^^ 

Alternatively, performing a good single-shot measurement in the sub-Zcno regime consists of three tasks: 

1. Localise the qubit eigenstates by turning on the qubit bias (e) and/or decreasing qubit tunnelling rate (A), so 

2. Increase the lead tunnelling rate (Tqq and Xoo) by lowering a tunnel barrier. 

3. Turn on the lead bias voltage (V). 

The rate at which 9 is varied effects a rotation on the qubit, so by selecting this rate appropriately, we can choose to 
measure in an arbitrary basis. 

XI. CONCLUSION 

In this paper we have used the quantum trajectories formalism to derive master equations for measurements 
of a charge qubit by an external point contact electrometer in the sub-Zeno regime. The master equation was 
derived without recourse to heuristic arguments, and resulted in the inclusion of inelastic processes in which energy is 
exchanged between the qubit and the detector. These inelastic processes have a profound effect on the conditional and 
unconditional dynamics of the system. Furthermore, within the sub-Zeno regime, our results are valid for arbitrary 
detector bias voltage. 

In the low-bias regime {eV < (f) the qubit always relaxes to its ground state, much like a qubit in equilibrium with 
a zero-temperature bath. In this case, relaxation to the ground state is due to the eventual spontaneous relaxation 
of the qubit accompanied by an excitation of a PC lead electron. The corresponding steady state power spectrum is 
flat. 

In the high-bias regime {eV > 0) the PC leads act like a zero-entropy heat bath at non-zero temperature and 
both inelastic relaxation and excitation processes take place, causing the steady state of the qubit to be a mixture of 
excited and ground states. These inelastic transitions are also reflected in the power spectrum, which exhibits a peak 
centered at zero frequency, corresponding to transitions of the qubit between the localized energy eigenstates. 

The sharp transition the dynamics at eV = (/) is also reflected in the steady state conductance of the PC. This ob- 
servation provides techniques for accurately determining the qubit Hamiltonian parameters directly from conductance 
measurements. 

Single shot measurement remains possible in the sub-Zeno regime, provided that the system Hamiltonian is modified 
such that the energy eigenstates become localized before the measurement takes place. This also gives the added 
freedom that we may choose an arbitrary basis in which to measure. 
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the E.U. NANOMAGIQC project (Contract no. IST-2001-33186). We thank H.-S. Goan, W. J. Munro, T. Spiller, 
G. J. Milburn, H-A. Engel, R. Aguado, A. Slmirman, H. M. Wiseman and D. Avcrin for useful conversations. In 
particular, H. M. Wiseman suggested Fig. 5(a). 

APPENDIX A: GENERALIZED POWER SPECTRUM OF GOAN AND MILBURN^ 

The power spectra calculated in^ assumed e = 0. In order to compare our power spectra with theirs for arbitrary 
parameter values, we derive the power spectra for their model for arbitrary e. 

It is laborious to compute Ggm(''") directly, however since we are only interested in the power spectrum, we can 
bypass the explicit solution for the correlation function in the time domain and calculate the power spectrum directly. 

We wish to compute the Fourier transform, F, of equation (A8) of^ (hereafter called (GMA8)). This may be done 
via a Laplace transform, L, using the relation for a symmetric function fit) 

F.[/(t)] = U^[/(t)] + L_,^[/(i)]. (Al) 

The Laplace transform of (GMA8) is 

U[Ggm(t)] = e/ioc + «i^oe(Tr{niU[e^^nipoo]} - Tiinip^Y / s). (A2) 
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The notation e^'^riiPoo is just shorthand for the solution to the unconditional master equation (GM5a) at time t 
subject to the initial condition p(0) = riipr^ = ni/2. Taking the Laplace transform of (GM5a) gives 

sLs[p{t)] - p{0) = -i[Hsys, LMt)]]+T^[Ti + {Tr - Ti)n{\LsW)]. (A3) 

Ls[p(<:)] may be found straightforwardly from this expression and then Ls[e'^'^ni/9oo] is obtained by substituting p(0) 
niPoa = 'T'i/2 into Ls[p(r)], though it is a somewhat cumbersome expression. Evaluating the traces in Eq. (A2) yields 

Ls[^gmIt;j 2 ^°'i°M(4s3 + 452^-2 +(/,2A'2 + s(4(/.2 + A'4)-(/)2;f 2 cos(20))' ^ > 

where we have used the same definition of X as found in Goan and Milburn^, which for comparison is given in our 
notation as X — v^/eV . The power spectrum is then given by 

Sgm{<^) = 2F^[GGM(i)] = 2(U[Ggm(t)] + L_,^[Ggm(t)]) (A5) 
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